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Weak-lensing 5-modes as a probe of the isotropy of the universe 
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We compute the angular power spectrum of the 5-modes of the weak-lensing shear in a spatially 
anisotropic spacetime. We find that there must also exist off-diagonal correlations between the E- 
modes, 5-modes, and convergence that allow one to reconstruct the eigendirections of expansion. 
Focusing on future surveys such as Euclid and SKA, we advertise that observations can constrain 
the geometrical shear in units of the Hubble rate at the percent level, or even better, offering a new 
and powerful method to probe our cosmological model. However, the power of this new technique 
still requires further investigations and a full analysis of signal-to-noise ratio. We hope that it can 
provide a new tool for precision cosmology. 

PACS numbers: 98.80.-k, 98.80.Es, 04.20.-q 


According to the standard lore p], in a homogeneous 
and isotropic background spacetime, weak-lensing by the 
large scale structure of the universe induces a shear field 
which, to leading order, only contains ill-modes. The 
level of 5-modes is used as an important sanity check 
during the data processing. On small scales, 5-modes 
arise from intrinsic alignments [2], Born correction, lens- 
lens coupling [3], and gravitational lensing due to the red- 
shift clustering of source galaxies |4|. On large angular 
scales in which the linear regime holds, it was demon¬ 
strated [5] that non-vanishing 5-modes would be a sig¬ 
nature of a deviation from the isotropy of the expansion, 
these modes being generated by the coupling of the back¬ 
ground Weyl tensor to the 5-modes. 

In this letter, we emphasize that, as soon as local 
isotropy does not hold at the background level, there ex¬ 
ist a series of weak-lensing observables that allow one 
to fully reconstruct the background shear and thus test 
isotropy. We also quantify their magnitude for typical 
surveys such as Euclid [7] and SKA [5j. As a conse¬ 
quence of the existence of 5-modes, it can be demon¬ 
strated that: (1) the angular correlation function of the 
5-modes, C BB , is non-vanishing [5j; (2) they also corre¬ 
late with both the 5-modes and the convergence, lead¬ 
ing to the off-diagonal cross-correlations (5^ m 5^ ±lm _ M ) 
and in which 5^ m and 5^ m are the com¬ 

ponents of the decomposition of the 5- and 5-modes of 
the cosmic shear in (spin-2) spherical harmonics, and 
are the components of the decomposition of the conver¬ 
gence in spherical harmonics [5]; (3) the deviation from 
isotropy also generates off-diagonal correlations between 
k and 5-modes, (E im E* t±2rn _ M ), {nf m nf ± * 2m _ M ), and 
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Our companion article [5] provides all the technical 
details of the theoretical computation of these correla¬ 
tors. In this letter, we estimate the information that 
can be extracted from weak-lensing by focusing on these 
correlations, and illustrate its power to constrain a late 
time anisotropy. 

We assume that the background spacetime is spatially 
flat and homogeneous, but enjoys an anisotropic expan¬ 
sion. It can be described by a Bianchi I universe with 
metric 

ds 2 = — df 2 + a(t) 2r )ij(t)<lx l &x : ’, (1) 

where a(t) is the volume averaged scale factor. The spa¬ 
tial metric 7^ is decomposed as 7 ij(t) = exp [2/3i(t)\dij 
with the constraint . =1 Pi = 0. The geometrical shear, 
not to be confused with the cosmic shear, is defined as 

a ij = ' (^) 

Its amplitude, a 2 = ay, a'i = Y^i=iPii characterizes 
the deviation from a Friedmann-Lemaitre spacetime. We 
also define H = a/a. 

At this stage, it is important to stress that, since a l:j 
is traceless, it has 5 independent components. In the 
limit in which a/H <C 1 (the relevant limit to con¬ 
strain small departures from isotropic expansion) each 
of the 5 correlators - {B lm E} ±l m _ M ), ( B im K * e±lm _ M ), 

(5^ m 5^ ±2m _ M ), (Kf m Kf±2m-M)’ an d {Bfrn K f±2m-M) ~ 

is of first order in <j/H and has five independent com¬ 
ponents (M = — 2... + 2) that allow one, in principle, 
to reconstruct the independent components of <Tjj. The 
angular power spectrum C BB , on the other hand, scales 
as (cr/H) 2 and, while it can point to a deviation from 
isotropy, it does not allow one to reconstruct the princi¬ 
pal axis of expansion. 
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Following our earlier works mm, we adopt an observer 
based point of view, i.e. we compute all observable quan¬ 
tities in terms of the direction of observation n 0 . The 
main steps of the computation are the resolution of the 
background geodesic esquation (which provides the local 
direction n(n 0 ,t ) on the lightcone and hence the defini¬ 
tion of the local Sachs basis), the resolution of the Sachs 
equation at the background level and at linear order in 
perturbation, and a multipole decomposition of all the 
quantities, a step more difficult than usual because of the 
fact that n ^ n Q . We then perform a small shear limit in 
which one can isolate the dominant, followed by the use 
of the Limber approximation (although not mandatory). 
This provides the expressions of the different correlators 


and L = £ + 1/2. P(fc) stands for the primordial 
power spectrum of the metric fluctuations, T v (x, rf) is 
the transfer function of the deflecting potential given, as 
usual, by the sum of the two Bardeen potentials. They 
are both evaluated on the past lightcone parametrized 
by the radial coordinate %; in the Limber approximation 
k = L/x■ Af(x.) is the source distribution and depends on 
the survey. a£ m (x) is the multipolar coefficient of the de¬ 
flection angle a(n 0 ,x) expanded in spherical harmonics. 
At the lowest order in a/H , only its l = 2 components 
are non-vanishing (see §VII.B.2 of [BJ for the expressions 
for a 2m ). 
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that take the form (see Eqs. (7.17-7.19) of Ref. [5]) 
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sF t l2 e 2 are explicit functions of the multipoles given in 
Appendix D.4 of Ref. (f> and defined in Ref. .TDj. The 
general form of the quantities is given by Eq. (7.14) 
of Ref. [6] and, in the Limber approximation, they reduce 
to 


FIG. 1: Time evolution of models A (red) and B (blue). The 
plot shows the contribution to the expansion of matter (solid 
line), dark energy (dashed line) and geometrical shear (dotted 
line, and magnified by a factor 100). 
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with (see Eqs. (7.20) and (7.21) of Ref. [B]) 


Vim = f -J^P (— ) &2m(x) 

Jo X Vx, 
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(5) 


We obtain 16] that for the R-modes 


While the previous off-diagonal correlators are the 
most direct consequence of a late time anisotropy, most 
experiments are designed to measure the angular power 
spectrum. 
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where 2 Fg 2 i+s is a function of t and the functions gf are expressed in terms of spherical Bessel functions and given 
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by Eq. (6.44) of Ref. [ 6 ]. 


To estimate these correlators, we need to go through 
the following steps. First, we need to solve the geodesic 
equation for the background spacetime in order to deter¬ 
mine n(n 0 , x) and the deflection angle. We then need to 
describe and solve the evolution of metric perturbations 
(in order to determine the transfer function T v of the 
lensing potential.) 

During inflation, the spacetime isotropizes, letting only 
tiny, if any, signatures on the cosmic microwave back¬ 
ground (CMB) pi], which has been constrained obser- 
vationally Hina. On the other hand, many models of 
the dark sector [ 12 ] have considered the possibility that 
dark energy enjoys an anisotropic stress. This is a generic 
prediction of bigravity [13] and backreaction |TB]. This 
has stimulated the investigation of methods to constrain 
a late time anisotropy, using e.g. the integrated Sachs- 
Wolfe effect [13], large scale structure and the Hubble 
diagram of supernovae in different fields mm- 

From a phenomenological point of view, one can con¬ 
sider a dark energy sector with an anisotropic stress. 
Its stress-energy tensor is then decomposed as T{? = 
(p + P)u^u v + P8^ + n(j where the anisotropic stress ten¬ 
sor n() is traceless (njj = 0 ), transverse (u^H^ = 0 ), and 
has 5 degrees of freedom encoded in its spatial part n*. 
It can be decomposed in terms of an anisotropic equation 

of state dill as P? = pdc (wSj + Awi'j. Here, w is 

the usual equation of state (we assume w = — 1 as for a 
cosmological constant) and we need to model H). 

The background equations then take the form 


3 H 2 

, . 1 2 
— K\Pm T Pde) 4 “ 2 ® ) 

(7) 

K)' 

= -3Ha) + kR) . 

( 8 ) 

Pm 

= -3 Hp m , 

(9) 

Pde 

= -cJijW*. 

( 10 ) 


The first equation is the Friedmann equation, the second 
is obtained from the traceless and transverse part of the 
Einstein equation and dictates the evolution of the shear. 
The last two equations are the continuity equations for 
the dark matter (P = H) = 0) and dark energy sectors. 

Simple models can be built by phenomenologically re¬ 
lating II) to the geometrical shear as If) = Act* = cr*/Kr n , 
where r n can be time-dependent. When it is constant, 
the shear grows exponentially as cr) = B* (^) e t/,T n. 
Since <J t j is small today, there is some fine tuning. We 
thus consider two classes of models defined by 

(A): Uj = pdeAwj ; (P) : A) = g{a)Aw). (11) 

This assumes that the anisotropic stress evolves with 
time while keeping its eigenvalues in a constant ratio. 
The function g(a) is arbitrary, and we note that when 


g{a) = 3H/H 0 , er* = Cj (yf ) 3 + so that at late 

time a 2 oc k 2 Aw 2 /Hq is constant. In the models (A), 
the dark energy triggers the anisotropic phase. It was 
argued [15] that next generation surveys will be capable 
of constraining anisotropies at the 5% level in terms of 
the anisotropic equation of state, a number to keep in 
mind for comparison with weak-lensing. Eq. ([ 8 ]) implies 
that 



ri ( d{a/a 0 ) 
j w H 


while Eq. (10 1 implies that pd e will decrease as 


Pde — PdeO exp 



Figure [T] depicts the evolution of the density parameters 
for a model of each class. 

To evaluate the angular power spectrum of the El¬ 
and B-modes, one needs to specify W(y). For this we 
consider the distributions of the future Euclid and SKA 
experiments. The normalised Euclid redshift distribu¬ 
tion [3 [ 2 D] is 


Af(z) = Az 2 exp 



( 12 ) 


with A = 5.792, /3 = 1.5 and z o = 0.64. For SKA, we 
use the SKA Simulated Skies simulations m of radio 
source population using the extragalactic radio contin¬ 
uum sources in the central 10 x 10 sq. degrees up to 
z = 20. The SKA normalized redshift distribution is [22] 


A f(z) = A 


(1 +z) r 


■ exp 


(a + bz ) 2 

1X+W\ ’ 


z<20 (13) 


with best fit parameters a = —1.806,6 = 0.388, to = 
2.482, n = 0.838 and A = 1.610, which give a description 
accurate to the percent level. 

Figure [2] depicts the two angular power spectra for 
these two surveys. In the linear regime the B-mode 
contribution is expected to vanish and the terms (. a Timi 
proportional to the off-diagonal correlators A A / ±1 , are 
shown on Fig. [3] However, the shear induces a B-mode 
spectrum whose amplitude is about ( a/H ) 2 lower than 
for the one for the B-mode in the most optimistic model 
(B). We can compare our results to the bounds set by 
the CFHTLS survey [22]. Unfortunately, CFHTLS cov¬ 
ers 4 fields of typical size 50 sq. deg. so that the largest 
scale with a sufficiently good signal-to-noise ratio is of 
the order of l ~ 2000, far beyond the linear regime. B- 
modes will be generated from the non-linear dynamics 
and it is safer to rely on the BB cross-correlation. To 
get a rough idea, though, we use the values at £ = 2000 , 
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FIG. 2: Angular power spectra of the E- and B-modes (resp. 
solid and dashed lines) for the Euclid (red lines) and SKA 
(blue lines) surveys for models A (long dashed) and B (short 
dashed). 


for which e-Cf E ~ 1CT 1 2 3 4 5 6 and £ 2 Cf B < 4 x 10" 7 . In¬ 
deed this estimation has to be taken with a grain of salt 
given the fact that i) there is a large scatter in Fig. 6 
and 7 of Kitching et al. (2014) and ii) these observations 
are not in the linear regime and there is no unambigu¬ 
ous way to scale these observations to lower £. We can 
thus set the bound ( a/H) 0 < 0.4 from CFHTLS. How¬ 
ever, Euclid shall probe scales up to 100 deg., deep in 
the linear regime, with a typical improvement of a factor 
50 [25] ■ This would translate to a sensitivity of order 
(a/H)o < 0.4/50 ~ 1% for the shear. This estimate 
indicates that weak-lensing could be a powerful tool to 
constrain a late time anisotropy. In the meantime, ex¬ 
periments such as DES will allow us to forecast more 
precisely the power of Euclid. It demonstrates that in 
principle one can reconstruct the principal axis of expan¬ 
sion from observations. We thus want to draw the at¬ 
tention on the importance of these estimators and their 
measurements. 

Conclusion: this letter emphasises the specific signa¬ 
tures of an anisotropic expansion on weak-lensing, as first 
pointed out in Ref. [5]. Following our formalism detailed 


in Ref. [B] (where all the technical details can be found), 
we have focused on two phenomenological anisotropic 
models and computed the angular power spectra of the 
E- and B-modes, as well as the five non-vanishing off- 
diagonal correlators. These are new observables that we 
think must be measured in future surveys. These mea¬ 
surements can be combined easily with the Hubble dia¬ 
gram since the Jacobi matrix can be determined analyt¬ 
ically at background level [3j. Let us emphasize that the 
off-diagonal correlations with the polarisation can also 
be applied to the analysis of the CMB, hence generalis¬ 
ing easily those built from the temperature alone [23]. 

Our analysis demonstrates that future surveys, and 
in particular Euclid, can set strong bounds on the 
anisotropy of the Hubble flow, typically at the level of 
(a/H)o < 1%, but one needs a detailed analysis of the 
signal-to-noise ratio to confirm this number, which for 
now has to be taken as an indication. This is a new 
and efficient method and the estimators built from the 
off-diagonal correlators can be used to reconstruct the 
proper axis of the expansion. 


'’VVw-Qi 



FIG. 3: The five correlators 70 y_i are, up to shape factors, 
proportional to CVim- We plot the latter for the models A 
(dashed) and B (solid) for the Euclid (red) and SKA (blue). 
For clarity we plot only the M = 0 component. 
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